The narrow escape problem revisited 
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The time needed for a particle to exit a confining domain through a small window, called the 
narrow escape time (NET), is a limiting factor of various processes, such as some biochemical 
reactions in cells. Obtaining an estimate of the mean NET for a given geometric environment is 
therefore a requisite step to quantify the reaction rate constant of such processes, which has raised a 
growing interest in the last few years. In this Letter, we determine explicitly the scaling dependence 
of the mean NET on both the volume of the confining domain and the starting point to aperture 
distance. We show that this analytical approach is applicable to a very wide range of stochastic 
processes, including anomalous diffusion or diffusion in the presence of an external force field, which 
cover situations of biological relevance. 

PACS numbers: 87.10.+e,05.40.Fb,05.40.Jc 
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The first-passage time (FPT), namely the time it takes 
a random walker to reach a given target site is known 
to be a key quantity to quantify the dynamics of vari- 
ous processes of practical interest [l|, 0,11]. Indeed, bio- 
chemical reactions 0, H, 0] , foraging strategies of an- 
imals [1, 0, [1(3] , the spread of sexually transmitted dis- 
eases in a human social network or of viruses through the 
world wide web [III ] are often controlled by first encounter 
events. 

Among first-passage processes, the case where the tar- 
get is a small window on the boundary of a confining do- 
main, defined as the narrow escape problem, has proved 
very recently to be of particular importance [T3 |. The 
narrow escape time (NET) gives the time needed for a 
random walker trapped in a confining domain with a sin- 
gle narrow opening to exit the domain for the first time 
(see fig HJ . The relevance of the NET is striking in cel- 
lular biology, since it gives for instance the time needed 
for a reactive particle released from a specific organelle 
to activate a given protein on the cell membrane [l3j |. 
Further examples are given by biochemical reactions in 
cellular microdomains, like dentritic spines, synapses or 
microvesicles to name a few [l2l [l3j . These submicrom- 
eter domains often contain a small amount of particles 
which must first exit the domain in order to fulfill their 
biological function. In these examples, the NET is there- 
fore a limiting quantity whose quantization is a first step 
in the modeling of the process. 

An important theoretical advance has been made re- 
cently by different groups [12, EH, HE] , which obtained the 
leading term of the mean NET in the limit of small aper- 
ture in the case of a brownian particle. However, these 
different approaches lose track of the dependence of the 
mean NET on the starting point. Obtaining such infor- 
mation is not only an important theoretical issue, but 
also a biologically relevant question. As a matter of fact, 
biomolecules like membrane signalling proteins or tran- 
scription factor proteins are generated at specific sites 
in the cell [l3j], whose localization plays an important 
role in the very function of the biomolecules, as under- 




FIG. 1: The narrow escape problem : the particle starts from 
rs and evolves in a domain Q with reflecting walls, except a 
small aperture S a of typical radius a centered at tt- 



lined recently in Q. In addition, the above mentioned 
techniques to estimate mean NETs have been limited so 
far to normal brownian diffusion, whereas many experi- 
mental studies have shown that cellular transport often 
departs from thermal diffusion due to the complexity of 
the cellular environment. In particular, crowding effects 
have proved to induce subdiffusive behavior in many sit- 
uations [IE, E3]> while the interaction of a tracer particle 
with molecular motors induces a biased motion [la ]. 

In this Letter we propose a general theory which (i) 
provides explicitly the scaling dependence of the mean 
NET on both the volume of the confining domain and 
the source-target distance, (ii) applies to a wide range of 
transport processes, including anomalous diffusion, (iii) 
encompasses the case of transport in the presence of a 
force field. Our formalism partially relics on the method 
recently proposed in [2], and considerably broadens its 
field of applications. 

We start by extending the theory developed in [2| to 
compute the mean FPT of a continuous random motion 
to a closed surface S a (jCT) of typical radius a containing 
tt, starting from a source point r$. We will next give 
explicit results when S a (rT) is a sphere, and show that 
this approach allows to solve the narrow escape prob- 



2 



lem. The volume delimited by So,(tt) will be denoted by 
P a (r^). We consider that the random walker evolves in 
a bounded domain £1 of volume V of the d-dimensional 
space M. d . Let P(r, t\r') be the propagator, i.e. the den- 
sity probability to be at r at time t, starting from r' at 
time 0, which satisfies the backward equation [18l |: 



d_ 

dt 



P(r,t\r') = A r ,P(r,t\r'), 



(1) 



where A r denotes the Laplace operator and the diffusion 
coefficient has been set to 1. Let F(r' ,t\r)dS(r') be the 
probability that the first-passage time at the infinitesi- 
mal surface dS(r') located at r', starting from r, is t. 
By partitioning over the first arrival time t' at a surface 
element dS(r') of the sphere S^ry), one obtains 

P(r T ,t\v s )= f dS(r) ( P(Y Tl t-t'\v)F{v,t'\Ys)dt'. 

Jr£S a (r T ) JO 

(2) 

We next assume a <C V 1 ' d , so that P(r T , t\r G S a (r T )) = 
P(r T ,t\S a (r T )) does not depend on r G S a (r T ). This 
condition will be fulfilled in the large V limit consid- 
ered in the following. Denoting p(s) = p(t)e~ st dt the 
Laplace transform of p(t), and setting 

F(S a (v T ),t\r s ) = fdS(T)F(T,t\v s ), (3) 



we get from equation @: 



P(r T , s\r s ) = P(r T , s\S a (r T ))F(S a (r T ), s\r s ). (4) 

We now denote by (T)(5 a (r T )|r s ) the mean FPT at 
S (r T ), and write finite P(r T ,*|r s ) = P Bta t(i*T)< Then 
one has 



F(S a (r T ),s\rs) = 1 - S (T)(S a (r T )\r s ) + o(s) 
P(r T , s\r s ) = P st at(r T )/s + H(r T \r s ) + o(l) 



(5) 



where 



poo 

H(r\r')= / (P(r,t\r')-P stat (r))dt. (6) 
Jo 

A similar expansion holds for P(rr, s|5 a (rr)). From Q, 
we then get 



<T)(S a (r T )|r s )P stat (r T ) = fl-(r r |S„(r r )) - H(t t \t s ), 



7) 



which is an extension of a similar form given in 0, [l 
We then consider the large volume limit of equation ([7|) 
and define the function a (rr|rs) as 



As can be checked directly from definition ([6]), (f> a (rT\rs) 
satisfies the following boundary value problem in the in- 
finite space: 



A r a (r T |r) = for r e K d \P Q (r T ) 
a (r T |r) = for r G 5 a (r T ) 



a„</» a (r T |r)d5(r) = 1. 

r£S a (r T ) 



(9) 



Equations (|8l9p completely define the large volume 
asymptotics of the FPT and constitute the central result 
of our method. Indeed, rephrased as above, the problem 
amounts to determining an electrostatic potential outside 
a conducting surface S a (rT) of charge unity, which has 
been extensively studied. When So,(tt) is a sphere the 
solution is straightforward and yields for the mean FPT: 



— \n(r/a) 
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lux , (T)/V = t\d,2) ( I 1 



2n d/2 



-,d-2 



„d-2 



for d 



for d > 3 



where r = r-r — r 
form found in 12 [t 



(10) 



;|. This result is compatible with the 
^211 using a different method. Besides 
the very useful analogy with potential theory, we will 
show in the following that the advantage of the formu- 
lation (|8l9p is threefold, (i) First, it can be adapted to 
other geometries and in particular to various examples 
of extended targets, such as an escape window in the 
domain boundary. It therefore extends the main result 
of [H obtained in discrete space for a point-like target, 
(ii) Second, as the derivation of equations (|8I9|) is inde- 
pendent of the operator A, it can be reproduced for any 
displacement operator C. In the general case, equations 
(|8l9p still hold, but with A to be substituted by the ad- 
joint operator C + , and d n 4> a to be substituted by the 
flux of <p a . (iii) Last, and following the previous remark, 
the formulation (|8l9p can be extended to the case of a 
random walker experiencing an external force field. 

We first show that equations (|8I9[) can be extended 
to the example of an escape window W a {vT) of small 
typical radius a centered at yt G 9f2, which is precisely 
the narrow escape problem. It will be useful to write 
Wai^r) = dfl fl B^{yt) where B^{tt) is a small volume 
of typical thickness e. We then set <S^(r-r) = QPidBKrr)- 
We now derive the mean NET through S*(tt), and we 
will use the fact that lim £ ^ S^t) — Wa( r T) to obtain 
the mean NET through W Q (rr). Following step by step 
the previous derivation, we obtain: 

(T)(^(r T )|r s )P stat (r T ) = ff(r T |^(r T )) - H(r T \r s ) 

We now take the infinite volume limit keeping tt — rs 
fixed, with the prescription that tends to the half space 
H.1 delimited by the hyperplane containing W a (rT), and 
define: 

hm (T)(^(r T )|r 5 )P stat (r T ) = ^(r T |r 5 ). (12) 

V— >oo 



1 / 1 ^ ^ T ^' 5 ' a ( rT ^ rs ) Pstat ( rT ) ~ M r r\rs)- (8) One can show that ^(r T |r) then satisfies: 
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A r #(r T |r) = for r e R d + \B e a (r T ) 
#(r T |r)=OforreS£(r r ) 
d n (f> e a (r T \r) = for r G dF d 



d n cj>l(r T \r)dS(r) 

&S>(r T ) 



-- 1. 



(13) 



This shows that 0*(rT , |r) is the electrostatic potential in 
a half space delimited by an isolating hyperplane contain- 
ing a conducting window <S^(r T ) of charge unity. Taking 
e to gives the mean NET through Wo,{yt)- In the case 
of a spherical window W a {vT), the solution of ([13]) can be 
exactly given. For d = 3, we obtain in oblate spheroidal 
coordinates 

lim (T)/V = — arctanf = —+o(-\ (14) 

where £ depends on cartesian coordinates according to 



z 2 x 2 + y 2 



= 1. 



a 2 ? 2 a 2 {e + l) 
For d = 2, we use elliptic coordinates and get 



(15) 



lim (T)/V = ± -r^eo - In(r/o) (16) 

V^oo 7T 7T 

where fj, depends on cartesian coordinates according to 



V 



a 2 cosh /i a 2 sinh /i 



(17) 



We stress that expressions (11411 6p of the mean NET are 
exact for any position of the sou rce p oint r <j . They there- 
fore extend the results of [H, EH Ell , which give the same 
small a limit, but where the dependence on the source 
position was not given. 

We now generalize equations (|8l9p to the case of a 
generic displacement operator C such that the station- 
ary distribution is uniform P s tat = 1/V, which actually 
underlies many models of transport in complex media 
[23| . We here assume that Sai^r) is a sphere, and fol- 
lowing 0, H3 , we further assume that the infinite space 
propagator Pq satisfies the standard scaling: 



Po(r,t|r')~* 



-d f /d m 



n 



(18) 



where the fractal dimension d f characterizes the accessi- 
ble volume V r ~ r d within a sphere of radius r, and the 

walk dimension d w characterizes the distance r ~ t 1 l dm 
covered by a random walker in a given time t. This for- 
malism in particular covers the case of a random walk on 
a random fractal like critical percolation clusters, which 
gives a representative example of subdiffusive behavior 
due to crowding effects 24| and could mimic in a first ap- 
proximation the cellular environment. Note that we here 
implicitly require that the trajectories and the medium 



have length scale invariant properties which ensure the 
existence of d w and df. Substituting the scaling |[TB]) in 
the definition ([5]), we obtain from ([5]) the large V equiv- 
alence: 

a ( a <u-d f _ r d w -d f ^ for dw < df 



lim (T)/V = { a\n(r/a) 

a ( r d-u,-d } _ a „ 



for d w = df 

' v 1 -<-dV-d' f „d w -d f ^ iordw>df 



Strikingly, the constant a does not depend on the con- 
fining domain but only on the scaling function II: 



IL(u 



-l/d v 



a = < 



o «' 
d w U{0) 

du 



d t /d u 



du 



df/d u 



n(o) -n(u- 1/d -)) 



for d u 
for d u 
for d n 



< df 
= df 
> df 



(20) 

Expressions (|19p therefore explicitly elucidate the depen- 
dence of the mean FPT on the geometrical parameters V 
and r. Note that for a generic target surface S a (rT), the 
r dependence in (flUj) still holds for r ^> a, but the a de- 
pendence is changed by a numerical factor. Indeed, the 
r dependence is fully determined by writing the conser- 
vation of the flux, which does not depend on the window 
shape at large r. As previously equations (|19p permit to 
obtain the mean NET: if we assume that the exit win- 
dow Sa(rx) is a half sphere of radius a, the mean NET 
will be exactly given by two times the mean FPT p9[) . 
For a generic window, and in particular in the case of a 
disk, the r dependence is unchanged for r 3> a, but the 
a dependence is modified by a numerical factor. 

Remarkably, equation (|19p highlights two regimes. 
When the exploration is not compact (du, < df), as in the 
case of a brownian particle in the 3-dimensional space, 
the dependence on the starting point disappears at large 
r. On the other hand, in the case of compact exploration 
(d w < df), as for 2-dimensional diffusion or subdiffusion 
on fractals, the mean NET diverges at large r and the 
starting point position is crucial. 

Last, we consider the case of a brownian particle in the 
presence of a force field F(r) = — V r $(r). In the context 
of biological cells, such force can be induced by the cou- 
pling of the particle to molecular motors which perform a 
directional motion along cytoskeletal filaments [l3j. We 
assume that the target is a sphere 5 a (rr), and that the 
force field is spherically symmetric and centered at r-r, 
situation which mimics asters of cytoskeletal filaments, 
which are ubiquitous in cells [HI . We set the gauge such 
that $(r) = for r G B a (r T ). Equations (EE]) then 
hold with A r to be replaced by the adjoint operator C + 
governing the evolution of the propagator [18| : 



C + = F(r)V r + A r 



(21) 



We then solve £ + (/> a (rT|r) = with the same boundary 
conditions as in ([5]), and write the stationary distribution 



(r) 



-*(r) 



/ n dfi(r') e-*C>- 



(22) 
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Using equation ||5J), we finally get the large V equivalence 
of the mean FPT: 

(23) 

As previously, the mean NET through a spherical win- 
dow is obtained as two times the mean FPT. One should 
note that the volume dependence entirely lies in the first 
integral factor of (f2"3")l . which is suitable for a quantitative 
analysis. On the other hand, the r dependence is fully 
contained in the second integral factor and will depend 
on the specific shape of $. The volume dependence of 
(|23p agrees with the one found in [25j, where however 
the r dependence was not given. 

We now give the explicit example of a divergence less 
force, which is a first approximation of the force experi- 
enced by a particle in an aster of cytoskeletal filaments, 
if one assumes that the force is proportional to the lo- 
cal concentration of filaments. This situation also de- 
scribes a brownian particle advected in an incompress- 
ible hydrodynamic flow. For d = 3 such force is given by 
F(r) — — 7/r 2 and the equivalent potential governing the 
dynamics can be taken as $(r) = j(l/r—l/a). Applying 
(23J>, we get 

ilmjT)/y = -L (er(V-Vr) _ ^ (24) 

for a generic domain shape. Note that the r dependence 
is modified by the hydrodynamic flow, while the V de- 



pendence at large V is unchanged for any flow intensity 
7. For d = 2 the force is given by F(r) = — 7/r and the 
equivalent potential can be taken as $(r) = 7ln(r/a). 
Applying (|23[) in the case of a domain whose boundary is 
parameterized by R(9) in polar coordinates, we get the 
large volume equivalence 

< T »~^(rs fl " (9, -Fi (r, - o,> ' 

(25) 

As opposed to the d = 3 case, the V dependence is now 
modified by the force. Interestingly we find a transition 
for 7 = 2. For 7 < 2, the mean FPT will scale like V 1 ^/ 2 
at large V, while the V dependence disappears for 7 > 2. 
In the case of a target centered in a spherical domain 
the results of [l[ can be straightforwardly recovered from 
(f2"3")l . and are compatible with (|24l25p . 

To conclude, we have proposed a general theory which 
provides explicitly the scaling dependence of the mean 
NET on both the volume of the confining domain and 
the source-aperture distance for a wide range of transport 
processes, including anomalous diffusion and transport in 
the presence of a force field, which are relevant to biologi- 
cal situations. In particular, we find that the dependence 
of the mean NET on the source aperture distance is cru- 
cial when the exploration is compact, as it is the case for 
the subdiffusive behavior in crowded environment which 
is observed in cells. Our model also permits to take into 
account the active transport due to molecular motors. 
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